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The behavior of the self diffusion constant of Langevin particles interacting via a pairwise in- 
teraction is considered. The diffusion constant is calculated approximately within a perturbation 
theory in the potential strength about the bare diffusion constant, ft is shown how this expansion 
leads to a systematic double expansion in the inverse temperature (3 and the particle density p. The 
one-loop diagrams in this expansion can be summed exactly and we show that this result is exact 
in the limit of small f3 and pf3 constant. The one-loop result can also be re-summed using a semi- 
phenomenological renormalization group method which has proved useful in the study of diffusion 
in random media. In certain cases the renormalization group calculation predicts the existence of 
a diverging relaxation time signalled by the vanishing of the diffusion constant - possible forms of 
divergence coming from this approximation are discussed. Finally, at a more quantitative level, the 
results are compared with numerical simulations, in two-dimensions, of particles interacting via a 
soft potential recently used to model the interaction between coiled polymers. 



I. INTRODUCTION 

Transport properties of deeply cooled liquids may 
change by several orders of magnitude when one re- 
duces the temperature [lj. fndeed, a liquid undergoing 
a deep quench below its melting transition may stay in 
a metastable supercooled state. In this state, the re- 
laxation time is much larger than the experimental time 
scale, and the system is always out of equilibrium. Ex- 
periments suggest that in "fragile" glass formers, there 
is a temperature To at which the relaxation time di- 
verges. Even if the mere existence of a divergence of 
the relaxation time r at finite temperature is still contro- 
versial, the experimental data can be often fitted by the 
so-called Vogel-Fulcher-Tammann law r ~ exp ^ T ^ Tg . 

Since such a glassy behavior has been observed in a wide 
range of materials, there have been huge efforts dedicated 
to computing the transport properties in supercooled liq- 
uids. Among them, the Mode Coupling Theory (MCT) 
0- C3- Efll appears to give results and predictions which fit 
remarkably well with data for many different systems [j| . 
However, despite its successes, the derivation of MCT ei- 
ther from the Mori-Zwanzig formalism £| or from fluctu- 
ating hydrodynamics remains quite obscure, and sys- 
tematic improvements seem difficult to implement and 
control. Hence, with the aim of understanding better 
the structure of the dynamics, alternative methods of 
computing transport properties of systems of many in- 
teracting particles such as supercooled liquids are worth 
examining. In this paper, we will develop a method to 
compute the long time self-diffusion constant which al- 
lows systematic double perturbative expansions, in the 
strength of the interaction and in the density of parti- 
cles. The approach is based on the Langevin dynam- 
ics for N particles with two-body interactions, though 



this approach may be generalised easily to three-body 
interactions. Such Markovian Langevin dynamics can 
be invoked in liquids over length and time scales where 
inertial effects become negligible. The Langevin equa- 
tion thus represents a course grained image of the system 
and the effective parameters and interactions used in the 
Langevin approach require microscopic derivation. The 
Langevin approach also naturally describes the dynamics 
of colloids in solution [?J, the Brownian noise is induced 
by the solvent and the effective interaction between parti- 
cles is composed of a pairwise direct interaction between 
the colloids plus additional hydrodynamic interactions 
induced via the solvent. 

We consider the interacting set of Langevin equa- 
tions for particles interacting via a pairwise potential 
V(Xi — Xj) depending only on the distance between the 
particles at temperature T in D dimensions: 

-A- = -X^dx-ViXi X,) + „?. (I) 

3 

The units of time are chosen such that the white noise 
field r]f has correlation function 

K{t)^{t'))=2n5 l3 5^5{t-t'), (2) 

where the angled brackets indicate averaging over the 
thermal noise. The term k is thus the bare diffusion con- 
stant of the particles in the absence of interactions. The 
fluctuation dissipation theorem or Einstein relation im- 
plies that X/k = l/T. This system of equations can be 
used to describe a colloidal system of interacting particles 
suspended in solution when hydrodynamic interactions 
are neglected. The neglecting of hydrodynamic interac- 
tions is justified where the direct two body interaction V 
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is of much longer range than the hydrodynamic interac- 
tions. A commonly cited example is charged colloids in a 
solution when the Debye length is very large with respect 
to the particle sizes. 

The effective macroscopic diffusion constant K e of par- 
ticle i is defined by 



lim (Xf(t)) = 2Dn e t. 



(3) 



It is the aim of this paper to develop a new technique 
for the calculation of K e . There are two basic routes to 
calculate K e . The first is based on the direct calculation 
of the diffusion constant USHEimmElH 
Pill IrH ITsl ligj . The second is based on the calculation 
of the modification of the response to a small external 
force on a given particle due the interaction with the 
other particles - the so called relaxation effect; the re- 
sulting value of n e is then determined from the Einstein 
or fluctuation dissipation relation |2(J, 1*2! lL I22L |23( . One 
approach to studying the dynamics of a tracer particle is 
to write an effective one particle Langevin equation with 
a non-Mar kovian memory kernel, derived via projection 
operator techniques llllllilllllllllllllllllilli. 
This kernel must then be computed by invoking approx- 
imation or closure schemes such as Mode Coupling-like 
approximations lid, cluster expansions or weak 
coupling expansions |3, El ■ Other approaches are 
based on closure schemes for the Smoluchowski equation 
0- IHl fill I21I |22| . which normally involve closing the 
hierarchy of equations for the joint probability density 
functions by replacing, for instance, the three body joint 
probability density function by its corresponding Kirk- 
wood superposition approximation. Both techniques can 
be handled to produce results which are exact to first 
order in the particle density o [HHl 111 ElE^ . 



II. DIAGRAMMATIC EXPANSION 

Here we use a technique based on a perturbative 
weak coupling expansion of the Smoluchowski or Fokker- 
Planck equation for the N particles in interaction. 
The form of this perturbation expansion is identical 
to that used to calculate the effective diffusion con- 
stant of a particle in a random potential. We denote 
by P(xi,X2-- - ,XAr,t) the probability density function 
for the particle displacements from their original posi- 
tions at t = at time t, that is to say the density of 
{X t (t) - Xi(0); 1 < i < N}. We take initial conditions 
Xi(0) = x- '' where the x- '' are independently and uni- 
formly distributed throughout the volume V of the sys- 
tem. In principal one could take other initial conditions 
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for the , notably one could take their equilibrium dis- 
tribution. However if the system is ergodic then the re- 
sulting behavior of the diffusion effective constant should 
be independent of this distribution. The disorder induced 
by the random initial conditions shows the link with dif- 
fusion in a random potential. The forward Fokker Planck 
or generalized Smoluchowski equation for P is 



dP 
~dt 



kV 2 P + AV • (PV0), 



(4) 



where V is the gradient operator on M. DN where D is the 
spatial dimension and N the number of particles. The 
potential 4> hi this formalism is given by 



(xi,X 2 • ■ • ,Xjv) = O (Xi+xJ O) ,X2+X^'' ■ • ■ ,XAr+X^ ; ), 

(5) 



,(0) 



where 



/> (xi,x 2 • • • ,Xiv) = ^2 V ^ ~ x j)- 



(0) 



From here on we shall denote by the vector x, without a 
particle index, the global position vector (xi,X2, • • -x^) 
in M. DN and by the vector k the corresponding Fourier 
vector (ki, k 2 , • • • kjv). If one defines 

P(k,s)= / dt [ dxexp(-si- ik- x)P(x,t), (7) 



it is straightforward to show that P(k, s) obeys 
1 



P(k, s) 



nk 2 + s 
A f dq 



nk 2 + s 



(27T) 



ND 



k-q^(q)P(k-q, S ) 



(8) 



We note that because of our choice of coordinates relative 
to the initial conditions, P(x, 0) = 5(x) and also 

0(q)-exp( iq .x( o ))^(q). (9) 

One has also that 

0o(q) = (2n)^ D Mci), 

with 



(10) 



Using these definitions one obtains the equation 



P(k,s) 



ik 2 



nk 2 



' / aTTd k ' iMq) exp(iq • x (0) )P(k - q, 



s). 



(12) 
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Eq. i|12|) can be solved iteratively leading to the expan- 
sion which is represented in Fig. QJ. We note here that 
this weak coupling expansion as it stands only makes 
sense for potentials which are bounded as V(x) is treated 
as a small perturbation. However the series can be re- 
summcd to obtain physical results for unbounded (such 
as hard core) potentials. In addition the behavior of soft 
potentials is of direct physical interest as they provide 
course grained descriptions of coiled polymers [25(, star 
polymers [2(| and micelles j2?J in solvents. 

Momentum is conserved at each vertex and the Feyn- 
man rules are 

• Each solid line (horizontal) corresponds to the bare 
propagator 



G (k,s) 



1 



ftk 2 + s 



(13) 



where k is the momentum in that line. 



• Each vertex (vertical wavy line) carries a factor 

-Ak-q VA)(q)exp( iq -x(°>), (14) 

where k is the in going momentum (from the left) 
and q is the momentum flowing into the vertical 
wavy line. 

• Each momenta q flowing into a wavy line is inte- 
grated over R DN with the measure dqJ(2w) D . 

We see that the ingoing and outgoing momentum on 
each of the generated diagrams is not conserved. This 
is because the spatial translational invariance of the sys- 
tem is not explicit. To study a spatially translational 
invariant system we average over the initial position 
x(°) throughout the volume V N . We shall take uni- 
form initial conditions although any initial conditions 
which have the property of spatial translational invari- 
ance should give the same asymptotic (late time) prop- 
erties for G(k, s) = (P(k, s)) . Here the angled brackets 
with the subscript indicate the average over the initial 
position vector x^ ) and is defined by 



(15) 



In the limit of large V this integration multiplies each di- 
agram by a factor of (2ir) ND S(J2 v q_ v )/V N where the q„ 
are all the momenta flowing into the upward wavy lines 
in each diagram at each vertex v. This momentum con- 
servation ensures that the momentum flowing into each 
diagram is the same as that flowing out and indicates the 
invariance by translation in space of the system averaged 
over its initial conditions. After taking this average over 
the diagrams in Fig. Q we obtain the diagrammatic ex- 
pansion for G(k, s). The Feynman rules are as before but 
with the following modifications 

• Each diagram carries an overall factor of 
(2ttY n ^ d /V n . 



• There are only n — 1 independent momenta for each 
diagram of n vertices by momentum conservation. 
Each of these momenta is integrated with the mea- 
sure d<\/{2ir) D as before. 

The form of the perturbation expansion shown in Fig. 
(|2J), shows that one can write G(k, s) 



G(k, s) 



1 



Kk 2 + s («k 2 + s) 



;D(k, s) (16) 



The term -D(k, s) can clearly be expressed in terms of 
one-particle-irreducible diagrams. Diagrams which are 
one-particle-reducible are those containing a propagator 
having a pure momentum k flowing through one of their 
propagators and can thus be factorised. The propaga- 
tor connecting two one-particle-irreductible diagrams is 
therefore the bare one l/(«;k 2 + s). If we denote the onc- 
particle- irreducible diagram contribution as S(k, s) then 

KK + S [KK Z + S) z [KK Z + S)° 

which sums to give 

1 



G(k, s ) 



«k 2 - E(k, s) + s 



(18) 



We now note that by the conservation of probability 
G(0, s) = 1/s and we thus expect that for small |k| 



G(k, s ) 



1 



Atk 2 - k 2 £(s) 



(19) 



where E(k, s) « k 2 E(s). In the limit of small |k| and s 
therefore 



G(k, s ) 



1 



(20) 



The effective diffusion constant for the particles is ex- 
tracted using the fact that 



2Dn e t= lim((X,(0-Xi^) 2 ), 



(21) 



where the angled brackets on the right hand side above 
indicate the average over the thermal noise and over the 
initial conditions. We also have that 

G(k, 8 ) = (exp (-i J2 • (Xi(t) - xf ^ ) (22) 
For small s in Laplace space we have 

/ dt exp(- S <)((X J (t)-x l (0) ) 2 ) 
Jo 

' 2Dk 
dt exp(-st)2Dn e t = — ^ (23) 



FIG. 1: Diagrammatic expansion for P(k,s) 
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FIG. 2: Diagrammatic expansion for G(k, s) (shown as a line with a blob) obtained after averaging over initial particle positions 
in the diagrammatic expansion shown in Fig. Q. 



and therefore, from Eq. (|22|l . for small s 



a— 1 * 



Now using Eq. H2U|I we obtain 

K e = k - E(0) 



(25) 



The calculation of n e can therefore be evaluated from 
G(k, 0). We may express the term E(0) as an expansion 
in the number of vertices v in the diagram, we write 



where 



E v = lim 0) 
|k|->0 k 2 



(26) 



(27) 



where S^(k) is the sum of one-particle-irreducible dia- 
grams with v vertices. From Fig. (J2J) we may write 

G(k ' 0) = i + 7^kW£^ (k ' 0) (28) 



where D v (k, 0) is the sum over all terms with v vertices 
and not just the one-particle-irreducible ones. The one- 
particle- irreducible components of D„(k, 0), E„(k, 0), 
must be extracted from -D„(k, 0). The terms Z?„(k, 0) 
have the behavior D v (\t,0) w F^k 2 for small |k|. It is 
straightforward to verify that the term F\ is zero. This 
means therefore that there is no contribution to O(A) 
in the asymptotic single particle diffusion constant K e . 
Examining Fig. J2J we find that 



£> 2 (k,0) 



{27r) n ( N-2) 

kV n 

k-q (k-q) 



(29) 



<iq 



(k - q) 2 



^o(q)^o(-q)- 



The potential term in the integrand may be expressed 
using Eq. as 

V-o(q)M-q) = E E v(<n)v{-<ik) 

i<j k<l 

S((h + qj)5(qk + (ll) X n 6 ^ n S ^ 

(30) 

One must keep in mind that only terms with nonzero q 
can contribute (one can imagine an additional term +.s 
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in the denominator for small s which is taken to zero at 
the end of the calculation) . The only terms in Eq. i|3U|) 
which have non-zero momentum are those where the pair 
= (k,l). In addition, the computation is simplified 
if one assumes k r = for r > 1 . Clearly the coefficient of 
in nk 2 — S(k, 0) is the self diffusion constant of particle 
1 which is also the self diffusion constant of any given 
particle. Hence for the purposes of the calculation of K e 
we can restrict our selves to the case k = (ki, 0, • • • 0). 
In this case, the only choices of the particle indices giving 
a nonzero diagram are i = k = 1, and j — I due to 
the scalar product k • q on the first vertex. There are 
thus N — 1 identical non-zero diagrams with two vertices. 
For notational simplicity in the following we will write 
k = ki eR D . All the choices are equivalent to choosing 
qi = q, q2 = q and q,. = for r > 2 where again 
q e M. D . In the resulting integral there are repeated 
delta functions for the terms r = s and also a double 



J 



5 (qi + <\j ) . We use the relation for q £ M. D 



(2tt 



(31) 



The result for Z?2(k) is thus 
£> 2 (k,0) 



\ 2 p f dq k ■ q (2q - k) • q ~ 
~J M » Wf (k-q) 2 +q 2 



The above diagram is also clearly one-particle- irreducible 
and hence F% = Ei thus giving 



1 \ 2 p f dq 



lX 2 p 



2 kD J rd {2-k) u 2 kD J rd 



The diagram giving Z?3(k, 0) has the value 



(33) 



( 2w )D(N-3) 



dq dp 



k ■ q (k q) ■ (q + p) (k + p) p 
(k-q)2(k + p)2 

I 



^o(q)^o(p)-0o(-q-p)- 



(34) 



To simplify the counting of diagrams with nonzero 
momentum let us consider the general expansion of 
V'o(qi) • • ' V'o(qn) in a diagram with n vertices. From 
Eq. (|lll) a given term on expanding the n-fold sum over 
pairs has the form 



A hjl (qi)---A injn (q n ). 



(35) 



Any diagram where the momentum q„ flowing into the 
vertex v is zero due to the presence of the scalar prod- 
uct with q„ at each vertex in the Feynman rules. The 
momentum flowing into the vertex v in the above de- 
composition over pairs is (0, 0, • ■ • q; • • • — q; ■ • • 0, 0) i.e. 
it has qi at the particle position i and — q^ at the particle 
position j in the total momentum vector. This is shown 
diagrammatically in Fig. 

Each horizontal line corresponds to a vertex and the 
horizontal coordinates are given by the points (i v , j v ). 
Each line must have a non zero momentum, thus q.„ 7^ 0. 
However the sum of the momenta down each column 
must also be zero. Hence each coordinate i v must ap- 
pear on at least two lines in order to give a diagram 
which is nonzero. In Fig. one example of the equiv- 
alent nonzero diagram contribution to E2 is shown. The 
3 types of diagrams which give nonzero contributions to 
F3 are shown in Fig. (J5J. The diagrams of type 1 have 
a multiplicity of iV — 1 and have the same particle pair 
on each line. For the diagrams of type 2 and 3 there are 
(N — 1)(N — 2) diagrams of the precise form shown in 
the diagram. This is because there are N — 1 choices for 
the first pair (1,2) and N — 2 choices of the pair (2, 3) 
(that is to say the particle number 3) shown on the sec- 
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FIG. 3: Diagrammatic representation of the vertex momen- 
tum present in one term of the pair development of a product 
of tpo's. 



ond line. The choice of the last (third) pair is not free 
by momentum conservation. Due to the presence of a 
delta function 5(J2 V Qi) f° r eacn column which is auto- 
matically satisfied when all the q^ are zero, each empty 
column carries a factor of V/(2ir) D . Hence a diagram 

with k empty columns has a factor ( (2n) D ) ■ 

grams of type 1, k = N — 2 and in diagrams of type 2 
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ko 



and 3, k = N - 3. 



q o o-q 

q 6 6-q 

The contribution of diagrams of type 1 to Ds(k, 0) is 
FIG. 4: Diagram in the pair development contributing to E2 thus 



n (l)f t nl X3 P f dq dp k ■ q [(2q - k) ■ (p + q)] [(k + 2p) ■ p] T >, xfr/ vfr , , > ^ 
^ (k '° ) = ~^L. (2^(2.)- [(k-q)' + q»] [(k + p)) 2+ p 2] nq)V(p)V(p + q), (36) 

which gives 



The contribution from diagrams of type 2 and 3 is These diagrams contributing to F3 are also one-particle- 

irreducible and hence to 0(A 3 ), we obtain 

1 



k 2iJfir Jrd \2tx) u zDk- 3 J R 2d (27T)-" (271")-^ p^q^ wn 6 J R n (2ir) 1J 

, (391 



Hence the ratio k & J k is expressed as a perturbation ex- 
pansion in y . Notice here that there are more lines than 
rows in the diagrams of the pair development, which 
means that a finite number of diagrams contribute to any 
given order of the ^-expansion, whereas an infinite num- 
ber of diagrams must be summed in order to compute any 
order in the p-expansion. A diagram with n lines and m 
non-empty columns is of order p" 1_1 (A/k)" = /9 m_1 /3", 
thus a systematic double expansion in p and j3 may be 
performed. We note that in Eq. (|39J) the leading order 
term in j3 which is of order pf3 2 recovers the weak cou- 
pling approximation. It is clear that the weak coupling 
approximation is not valid at low densities if the tem- 
perature is too low. The calculation of n e to first order 
in p involves summing all the diagrams which have just 
two columns occupied. These diagrams, which are all 
one-particle-irreducible, can be re-summed via an inte- 
gral equation [24| and one can show that the resulting 
expression for n e is the same as that given by the relax- 
ation method applied to the effective two-body problem 
as expounded in p^ . 



one particle index in common with the v — v' remaining 
vertices then the momentum flowing between the vertex 
v and v + 1 is zero (or k if the column in common is the 
first one) and thus the diagram is zero (or one particle re- 
ducible). In other words, all diagrams with disconnected 
loops such as the one shown in Fig. © have a zero value 
or are one-particle-reducible and thus give no contribu- 
tion to K e . 

As an example, we consider potentials of the form 

l ' (r) =(if)T e>q> (-€)- (40) 

which has been recently proposed to model the effective 
interaction between coiled polymers |25| at weak dilution. 
The Fourier transform of V(r) is then given by 

y(q)=er D exp^-^. (41) 

With this interaction potential and setting e = r = 1 we 
find from Eq. 1(55)1 



Consider a diagram with v vertices and its ipo pair ex- 
pansion. If the first v' vertices do not contain more than 



r q u pqP ± i 



(42) 
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FIG. 5: Diagrams in the pair development contributing to E3 (left: type one, center: type 2, right: type 3) 
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FIG. 6: Example of a diagram with disconnected loops in the 
pair development. Here the q-loop is disconnected from the 
(p,r)-loop. 



where the first non zero a BO 's are 



1 



012 

ai3 

«23 



1 / 1 K(D) + DK(D + 2) 
3 



(27T) 



D 



(43) 



K(D) = 



4 j D(6tt) i V 2 ' 
ax ■ 



l + x 2 



Clearly this y -expansion is valid only at very high tem- 
perature. It can be improved in several manners. One 
could compute higher orders in this expansion, which 
should give better agreement with the simulations, but 
however the expansion will inevitably break down at low 
T. Alternatively, one could try to sum infinite sub-series 
in order to build approximate non-perturbative schemes. 
As mentioned above one approach is to try and re-sum 
the diagrams to obtain results exact for all (3 to order 
p. Here we shall concentrate on another re-summation 
involving only one-loop diagrams. 



III. ONE-LOOP ANALYSIS 
A. Simple one-loop contribution 

Here we will focus on the class of one-loop diagrams. 
These diagrams are those which involve only one mo- 
mentum integral. Is this case, there are two dots on 
each occupied line or row and these diagrams are all one 
particle-irreducible. In addition from the discussion in 
the previous section, these diagrams are the dominant 
ones in the limit where pj3 = c (with c a constant) and 
p — ► 00, or equivalently where p/3 = c and (3 — > 0. The di- 
mensionless form of c would in fact be c' = pr® f3e, where 
ro is the characteristic length scale of the potential and 
e its energy. As before, we write k = ki e R D . The 
interesting point about this limit is that the (high tem- 
perature) statics of the model can also be evaluated [2^] , 
only chain diagrams in the virial expansion are retained 
in this limit. In electrolytic systems the Debye-Hiickcl 
approximation is recovered on the retention of only chain 
diagrams [2{j and hence it is interesting that one can have 
a theory of n e for electrolyte systems which is compatible 
with the Debye-Hiickel approximation which has proved 
so useful in the study of their static properties. 

The typical one-loop diagrams can be reduced to the 
staircase like diagrams shown in Fig. (JJJ by relabelling 
the particles. We note here that as we are summing an 
infinite number of diagrams one must be careful not to 
include the same particle twice in the same diagram as 
this will make the diagram two-loop. Hence for a finite 
system one cannot have a one loop diagram with more 
than N vertices as it will include at least one particle 
(other that the tracer particle (1)) at least three times, 
meaning that the diagram can have at least two two in- 
dependent momenta flowing through it. In this case the 
counting of the one-loop diagrams with n vertices where 
n is of order N will be different. We have however taken 
the limit TV — > 00 already and hence this does not cause 
us any problems in the region where the power series 
giving the one- loop result is convergent. 

In constructing a one-loop diagram at any level one 
must insert two points one of which must coincide with 
one of the two unpaired (in the vertical direction) points 
above. The leftmost point of the new pair can either 
be paired with the rightmost of the two unpaired points 
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FIG. 7: Example of a staircase diagram contributing to the one-loop expansion. 



above or it can be paired with the leftmost of the two 
unpaired points above - this gives a crossing. Clearly 
the first crossing must be back to the first column. Each 
diagram thus has p crossings where p G [l,n — 1]. If 
there is only one crossing then it must occur at level n 
corresponding to a staircase diagram with one crossing 
at the last line. Up to particle relabelling there is only 
one such diagram and its contribution is 



the second occurrence of the particle 1) and are given by 



b 2 {ni,n) 



2 n ~ 2 D k™- 1 



rfq 



■V(q) n . (45) 



Hence, taking into account all the multiplicities, the 
contribution of all the diagrams with n > 2 lines is 



2n- 3 p n_1 (-A) T 
2 n ~ 1 D K™" 1 



dq 

d (2^ 



V(q) n , (44) 



If n > 2 and p > 1 then the first crossing must occur at 
the level n\ S [2, n — 1]. At the level n\ + 1 the next two 
points can either cross or not cross just till the level n 
where the diagram most close. There are thus 2™~ 1 ~™ 1 
topologically distinct diagrams with the same first cross- 
ing point at n\. It is easy to see that these diagrams 
depend only on the position ri\ of the first crossing (i.e. 



J 



b{n) = h(n) + 2 n - n i- 1 b 2 (n 1 ,n) 

ni=2 

_ 1 - 2 1 -" p n - l {-\) n f dq 



(46) 



D 



n — 1 



d (2tt) j 



v( q y 



At n — 2 the above agrees the result of Eq. (|33|l . and 
hence the above formula is valid for all n>2. Performing 
the sum we find the one-loop contributions to be 



_ (one-loop) 



- K =-JZb(n) = -g- [ 7^(q) 2 ( ^ 



1 1 



V(q) 2i + £Ay (q); 



(47) 



The above may be conveniently rewritten in terms of c = p\/n 



K (one-loop) _ K 



l + ^( fl (c)-2 ff (£))' 



(48) 
(49) 



9 



where 



9(c) 



dq cV(q) 
o (2n) D l + cV(q) 



(50) 



From the previous discussions the corrections to Eq. I|49|) 
of the form 



Ke = K (™ c - lo °p) + /5 2 S2 (c)+/3 3 S3 (c)-. 



(51) 



In the particular case of the potential given by Eq. 
1 141)1 . we find that in D dimensions 



9(c) 



1 



(2^)-r(f) 



du 



cexp(-a)«2 1 
1 + cexp(— u) 



(52) 



In the case D = 2, we obtain 



(one-loop) 



Anp 



ln(l + p(3) - 2 In 



p/3 

1 + T 



However, this expression becomes negative when (3 is 
large, although this is outside the range where the ap- 
proximation involved here is valid. We also see that n e is 
a non-monotonic function of p, having a minimum value 
at some value p c but increasing up to k again on taking 
p very large - this is a consequence of the use of a soft 
potential, at high densities the particles all overlap but 
the energy change in moving to overlap with one particle 
rather than another is zero. The effective potential seen 
by the particles is almost flat and hence there is only a 
small effect on the particle diffusion. 



the renormalization-group method. This approach has 
been proven to be very accurate in the calculation of the 
effective diffusivity in random media and removes the fic- 
titious vanishing of n e predicted by simple perturbation 
theory. Indeed, in dimension three the one-loop rcnor- 
malization group analysis provides a very good quantita- 
tive approximation of the effective diffusivity of a p arti- 
cle in a Gaussian random potential |3Ct I'ilL l32l l8.'l 
and the exact result in dimensions one and two [35| • We 
therefore apply the same technique to the problem of 
interacting particles studied here. The potential tf> is de- 
composed into a long and short scale components 



<Mx) = 



dq 



q \<A (2tt) 
f/q 



|?|>A ( 27r ) 



\ D v(q) exp(iq • x) 
N£ >^(q)exp(iq-x) (54) 



where A is a running cut-off. One then integrates out the 
high momentum component 0> perturbatively to find an 
effective theory on length scales greater than 1/A. One 
then makes a self similarity ansatz which means only 
keeping interactions that were in the original problem 
and hence the only parameters which change are k and A 
(as well as the potential which just has the higher Fourier 
modes removed). One therefore has a running diffusion 
constant «(A) and a running coupling to the gradient 
field A(A). The effective diffusion constant is then given 
by K e — k(0). The flows of these couplings can be com- 
puted from the one-loop diagrams. In addition it can be 
shown by general arguments |35j that 



B. One-loop renormalization-group 

We have seen in the precedent section that the one- 
loop calculation of K e predicts that the diffusion constant 
can vanish at finite temperature. The same calculation 
predicts a vanishing diffusion constant for the problem 
of diffusion in a Gaussian random potential with short 
range correlations. This transition can be shown to be 
absent in finite dimensions via exact results in one and 
two dimensions, numerical simulations and general ar- 
guments based on the fact that the system has a finite 
correlation length. In order to go beyond the simple one- 
loop contribution, one can look at how the effect of the 
interactions on the diffusion constant propagates from 
short to large length scales. This can be achieved by 



«(A) K 
A (A) A 



(55) 



that is to say that the Einstein relation or Fluctuation 
dissipation relation is satisfied by the renormalized the- 
ory at each step of the renormalization. This rcnormal- 
ization is only valid for the low-momentum component 
of the remaining drift V</>< and a possible improvement 
to the calculation here would be to functionally renor- 
malize the field </><, which amounts to introducing new 
interactions generated by the renormalization procedure 
(this approach has been applied to diffusion in an incom- 
pressible quenched Gaussian velocity field j^ij) . Using 
the one-loop diagrams calculated in the previous section 
we find the flow equation for k(A) is 



dK (o noloop) (A) _ pX{A) 2 .. . ■ x x x , 

dA D K (A)(2n)° bD - lA V[A) [ 1 + £ MgLv {A) 2 1 + f ( A) 1 ' 
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where Sd-i is the area of the unit sphere of M. D . Insert- 
ing into Eq. (|56*|) one gets 



^(one-loop) ( A ) 



r/A 



PA(A) 
L>T(2vr) J 



1 + f V(A) 2 i + &V(A) 



(57) 



Integrating out A down from oo to and using the ini- 
tial conditions (k(oo), A(oo)) = (k, A), we find the final 
expression for the effective diffusion constant is simply 



K 



(one-loop RG) 



«exp 



K, 



(one-loop) 



(58) 



Now looking at K <° nc " loo P) gj ven j n Eq (|47[1 . the value 
of the self-diffusion constant given by the self-similarity 
ansatz leads to several comments: 

• The temperature-density dependence of n e / n is of 
the form h(^)/T. Hence the phase diagram in the 
(T, p) plane obtained from it consists of regions sep- 
arated by straight lines crossing at the point (0, 0). 
However, the Langevin approach studied here is in 
most physical situations a coarse-grained approach, 
in which case the potential V has to be replaced by 
an effective potential V e s (x) = V e g(x, c, p) which 
depends on the temperature and possibly the den- 
sity. In this case, the phase diagram will be more 
complicated. 

• If the Fourier transform V(q) of the potential 
has an absolute minimum of negative value at 
some non-zero value q* , then the diffusion con- 
stant goes to zero when the temperature reaches 
T* = —p(V(q*)^ 1 from above. Keeping the den- 
sity fixed lowering T amounts to increasing c. The 
integral 



1= dq 



1 



=V(q) 



(59) 



diverges at c = c* = p/T* , and near c* the denom- 
inator of the integrand behaves as 



(1 + (c* -Sc)(V(q*) + V qq (q*)Sq 2 /2)) 
= (-6cV(q*) + V qq (q*)Sq 2 /2), 



(60) 



where we have assumed that V is twice differen- 
tiable about q* , i.e. the minimum is not a cusp. In 
this case the relaxation time diverges at T* as 



exp 



A 



(61) 



However one could conceive that V is an effective 
potential coming from a course grained approach 



to a microscopic model and that it may thus have 
a dependence on p or T. Alternatively one might 
argue that diagrams of lower order in p dress the 
effective interaction leading to and additional p or 
T dependence in V . In this scenario, if we keep 
p fixed, the denominator of the integrand now be- 
haves as 

(1 + (c*-5c)(V(q*,c*)-V, c (q*,c*)5c+V, cc (q*,c*)Sc 2 /2 
+ V qq (q*,c*)Sq 2 /2 + V qc (q*,c*)5c5q)) 
« {5c/c* - c*V c {q\c*)5c-V qc {q\c*)8c5q)) (62) 

By definition we have that V c (q*,c*) < as the 
minimum of V(q,c) arrives at — 1/c from above. 
Hence is this case the relaxation time diverges as 



exp 



.4 



(63) 



which is the Vogel-Fulcher-Tammann law. From 
the above analysis we sec that if the dependence of 
V(q,c) on c is weak near (g*,c*) then one would 
expect to see a crossover between the behavior of r 
between Eq. (JHTJ for \T — T*\ < 1, to the behavior 
Eq. (E31> when \T - T*\ < 1. If q* = then 
the behavior is different and depends on the spatial 
dimension D. Finally as mentioned in the previous 
section on the one-loop expansion, precisely at T = 
T* one expects that K e depends explicitly on N. 

• In the corresponding static approximation |28| it 
was shown that the structure function is given by 



S(q) 



1 



l + cV(q)' 



(64) 



although the authors of |23 express S in terms of 
the Mayer function of the potential V, the result 
is the same at the order of accuracy of the cal- 
culation. We see that the condition above for a 
dynamic transition, defined by the vanishing of re e , 
coincides in the same approximation scheme by the 
appearance of a diverging correlation length and a 
second order phase transition. However as pointed 
out in |28j this apparent second order transition 
may be preceded by a first order freezing transi- 
tion when the free energy of the ordered crystal 
phase is lower than that of the free energy of the 
liquid phase described by the terms in the chain re- 
summed virial expansion. The analysis above thus 
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may apply to a supercooled liquid and the tempera- 
ture T* is the limit of the thermodynamic stability 
of the liquid phase. There has been some experi- 
mental evidence of a diverging correlation length at 
the Vogel-Fulcher-Tammann temperature T* from 
the study of the dielectric susceptibility of super- 
cooled liquids §57], though these results seem at 
odds with earlier numerical studies [38|, y!j study- 
ing this question. 

IV. COMPARISON WITH MONTE CARLO 
SIMULATIONS 



the one indicated as type 1 in Fig. (0 ■ The effect of the 
addition of this diagram to the one-loop RG calculation 
is plotted in Fig. JHJ as well. The other diagrams of 
the two- loop expansion are at least of order /3 4 . Hence, 
the good agreement observed indicates that the two-loop 
calculation which includes additional interactions will im- 
prove the approximation to the diffusion constant. We 
note that the value of the density used here (p = 0.5) is 
not very small, and the reasonably good agreement ob- 
tained is a hint that the RG calculation catches the p 
dependence of the diffusion constant in the case of soft 
or bounded potentials. 



In order to test the accuracy of the different schemes 
explained above, we have carried out Monte Carlo sim- 
ulations of particles interacting via the potential l|40[) in 
two dimensions. Here we have fixed the density at the 
value p = 0.5, and the number of particles at TV = 10000. 
We have evaluated the macroscopic diffusion constant by 
using Eq. The result is shown on Fig. JSJ and com- 
pared to the evaluation from Eq. <|42|l and <|431 for values 
of the interaction in the range [0, 7.5]. 



1.1 
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' • • • • 
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FIG. 8: Microscopic diffusion constant as a function of 
the interaction from different approximation schemes (a,b,c), 
compared to the one measured in Monte Carlo simula- 
tions: ^-expansion from Eq. ll 12t and l|4.SI (a), one-loop 
renormalization-group (b), one- loop renormalization-group + 
first two- loops diagram (c) and Monte Carlo simulations (d). 

From Fig. (jHJ we see that the /3-expansion is valid only 
at very small /3, however the one-loop RG result is much 
nearer the simulated values of the diffusion constant over 
quite a broad range. We expect that a RG analysis in- 
cluding more interactions, for example including func- 
tional renormalization of the interaction V, would pro- 
vide better agreement. The two- loop diagrams include 



V. CONCLUSION 

We have developed a perturbative expansion of the 
self-diffusion constant of Langevin interacting particles. 
The expansion is a weak coupling expansion which is 
suitable for soft or bounded interaction potentials. The 
perturbation expansion can be seen to be a double ex- 
pansion in [3 and p, and we have shown how it can be 
dealt with diagrammatically. As an example of a partial 
re-summation of this expansion, a one-loop renormaliza- 
tion group analysis has been carried out and tested on a 
simple form of the interaction, and compared to numer- 
ical simulations. It was found that this calculation gives 
considerably better agreement than the straight expan- 
sion to 0(/3 3 ). In addition, in some cases the one-loop 
RG calculation predicts a divergence of the relaxation 
time (or vanishing of the self diffusion constant) at some 
positive temperature. However we have seen that when 
such a dynamic transition occurs it is accompanied by 
a diverging correlation length in the statics when one 
uses equivalent approximations in the statics and dynam- 
ics. The possible forms of the divergence of the relax- 
ation time have been discussed and it has been argued 
that the Vogel-Fulcher-Tammann law emerges under rel- 
atively weak additional assumptions to the basic calcu- 
lation carried out here. The two-loop calculation should 
provide a better quantitative approximation for the dif- 
fusion constant, but it would also indicate the form or 
robustness of the divergence of the relaxation time as 
higher- loop contributions are taken into account, this 
calculation is in progress. Various other re-summation 
methods may be developed on the systematic perturba- 
tion expansion expounded here, as we have mentioned 
earlier the series may be re-summed at order p to re- 
cover existing low density results, one may use self con- 
sistent perturbation theory and also explore renormal- 
ization group schemes based of calculating the effect of 
integrating out the effect of a small density of other par- 
ticles rather than the Fourier modes of the interaction 
potential. 
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